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ABSTRACT 


Transition  curves  between  flat  faces  and  parallel  middle  bodies  are  developed 
as  families  of  special  polynomials  termed  ’’cubic”  polynomials.  The  curves  start 
and  end  with  zero  curvature  to  provide  no  discontinuities  in  curvature  at  the  junc¬ 
tion  with  the  flat  faces  and  the  parallel  middle  bodies.  Each  ’’cubic”  polynomial  is 
a  linear  combination  of  independent  polynomials  controlled  by  adjustable  para¬ 
meters.  Permissible  ranges  of  the  adjustable  parameters  are  examined  with  re¬ 
spect  to  selected  geometrical  constraints  such  as  inflection  points. 


ADMINISTRATIVE  INFORMATION 

The  work  described  in  this  report  was  sponsored  by  the  Naval  Ordnance  Sys¬ 
tems  Command  (Code  035B)  and  was  funded  under  UR  109-01-03. 

INTRODUCTION 

The  geometrical  characteristics  of  streamlined  bodies  of  revolution  are 
developed  in  Reference  l#by  the  method  of  independent  polynomials  and  adjustable 
parameters.  For  the  case  of  flat-faced  forebodies,  a  ’’cubic”  polynomial  is  intro¬ 
duced  for  the  transition  curve  from  the  flat  face  to  the  section  with  maximum 
diameter.  The  ’’cubic”  polynomial  has  the  desirable  property  of  providing  infinite 
slope  and  zero  curvature  at  the  junction  with  the  flat  face.  The ’’cubic”  polynomial 
of  Reference  1  has  one  adjustable  parameter  kj  ,  the  curvature  at  the  junction  with 
the  section  of  maximum  diameter. 

It  is  now  desired  to  develop  ’’cubic”  polynomials  that  fair  into  parallel  middle 
bodies,  that  is,  having  zero  curvature  at  the  junction  with  the  parallel  middle  body. 
For  the  ’’cubic”  polynomial  of  Reference  1,  kt  is  zero  to  satisfy  the  noncurvature 
of  the  parallel  middle  body  and  the  resulting  polynomial  has  no  adjustable  para¬ 
meter.  It  is  now  proposed  to  develop  ’’cubic”  polynomials  to  fit  parallel  middle 
bodies  so  that  the  polynomials  have  adjustable  parameters.  Two  adjustable  para¬ 
meters  are  considered: 

♦References  are  listed  on  page  17. 
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1.  ,  the  rate  of  change  of  curvature  with  arc  length  at  the  Junction  with  the 
flat  face. 

2.  ,  the  rate  of  change  of  curvature  with  arc  length  at  the  junction  with  the 
parallel  middle  body. 

The  polynomials  are  subjected  to  geometrical  considerations  such  as  those 
concerning  inflection  points  to  provide  suitable  ranges  of  values  for  the  adjustable 
parameters . 

The  pressure  distribution  on  a  body  in  a  flow  is  governed  by  the  curvature  of 
the  body  shape,  among  other  factors.  Discontinuities  or  sharp  changes  in  curva¬ 
ture  are  in  general  to  be  avoided  since  they  lead  to  pressure  fluctuations  which  may 
have  undesirable  results  such  as  separation  or  cavitation.  Such  changes  in  curva¬ 
ture  may  occur  at  junctions  of  bodies  or  at  inflection  points.  Hence  curvatures  at 
junctions  should  be  matched  and  inflection  points  avoided  in  most  cases. 

An  earlier  study  of  elliptic  curves  bridging  the  flat  face  and  the  parallel  middle 
body  is  given  in  Reference  2. 

TWO-PARAMETER  "CUBIC"  POLYNOMIALS 

The  "cubic"  polynomial  for  the  transition  curves  of  forebodies  between  flat- 
faced  noses  and  cylindrical  parallel  middle  bodies  requires  a  zero  curvature  at  the 
junctions.  For  the  "cubic"  polynomial  of  Reference  1,  =  0  and 

y3-  Q[x)  =  1  +  (x-1)3  (1) 

This  may  be  termed  the  zero-parameter  "cubic"  polynomial.  The  normalized 
coordinates  x  and  y  are  given  by  (see  Figure  1), 

x  =  X/X 

n 

y=  <2Y-Df)/(D~Df) 


2 


where 

X  is  axial  coordinate  starting  from  the  flat  face, 

Y  is  radial  distance  from  centerline, 

Xn  is  axial  length  of  forebody, 

D  is  diameter  of  parallel  middle  body,  and 
Df  is  diameter  of  flat  face. 

To  provide  a  more  general  family  of  bodies,  two  additional  adjustable  para¬ 
meters  are  introduced.  Since  curvature  k  represents  the  rate  of  change  of  slope 
with  arc  length,  the  rate  of  change  of  curvature  with  arc  length  £  ought  to  be  a  use¬ 
ful  parameter.  The  two  adjustable  parameters  are  then  and  ,  £  at  x  =  0  and 
at  x  =  1,  respectively. 

dk 

The  rate  of  change  of  curvature  with  arc  length  is  given  by  —  . 


"CUBIC"  POLYNOMIAL 


NORMALIZED  *.  ?Y  °t 

COORDINATES:  '  D  -  Df 

Figure  1  —  Basic  Geometry  of  Flat-Faced  Nose 
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Curvature  k  may  be  written  for  convenience  in  the  form 


which  yields 


At  x  =  0,  we  require 


—  =  d!* 

dy  dy* 


0 


and  then 


(T  = 


0 


Curvature  k  may  be  written  now  in  the  alternative  form 

-  (£)  [>*  te/]"* 


(where  x  is  taken  as  the 


independent  variable  since 


dx 


0  at  x  = 


9  dk  _  dk  dx 
ds  dx  ds 


(2) 


(4) 

(5) 

1) ,  which  yields 

>♦(£)]”  <6’ 
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At  x  =  lf  we  require 


and  then 


To  obtain  the  two- parameter  "cubic"  polynomials ,  the  method  of  Reference  1 
is  followed. 

There  are  two  conditions,  and  at 2,  for  the  two  adjustable  parameters: 


The  boundary  conditions  3 


j 


Pi: 

x  *  0, 

y  =  0 

x  =  1, 

y  =  i 

Bs: 

x  =  1, 

f- 

B4: 

X  =  1, 

S-. 

are 


Since  there  are  six  conditions  in  all,  the  required  polynomial  is  of  the  fifth 
degree . 
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The  and  0  ^  are  substituted  Into  the  "cubic"  polynomial 


y3  =  a  +  a4x  +  a2X2  +  a3x3  +  a4x4  +  a^5 


(8) 


to  give  for  each  condition 

6 

ai 

k 

o 

cr  2:  2aa  +  8a4  +  20a5  =  kt 
®ls  ao=0 

02*  aQ  +  aj  +  a2  +  a3  +  a4  +  a5  =  1 

03:  at  +  2a2  +  3a3  +  4a4  +  5a5  =  0 

0  4;  2a2  +  6&J  +  12a4  +  20a3  ~  0 

The  form  of  the  "cubic"  polynomial  is  then 


y3  =  “■  Kq  [x]  +  k,  K,  [x]  +  Q  [x]  (9) 

k 

o 

where  Kq,  Kj  and  Q  are  fifth-degree  polynomials  which  are  determined  as  follows 

The  preceding  relations  for  at  ^  and  0^  correspond  to  conditions  on  y3,  the 
"cubic"  polynomial,  instead  of  y,  as  follows: 


.  dy3  frti  6 


at:  [l)  =  6as  +  24a«  +  60a5 


3k, 


B,:  y3  [0]  =  aQ  =  0 
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02 :  y3  Cl]  =  a0  +  at  +  a2  +  a3  +  a4  +  a5  =  1 
Bj*  3x^  til =  ■*■  2&j  +  3ftj  +  4a*  +  5flj  =  0 

j2  3 

0^:  — [l]  =  2a2  +  6a2  +  12aj  +  20a§  =  0 

In  terms  of  the  polynomials  Lx]  ,  Kt[x3and  Q[x]the  preceding  relations 
represent 

crt:  K  ’[0]  =  6,  Kt’  Co]  =  Q’  Co]  =  0 

o 

or*:  K,’"  [l]=  3,  K  [l]=  Q"'  Cl]=  0 

P«:  K  10]=  K«  [0]=  Q  [0]  =  0 

o 

Pi:  KQCl]=  KjC  1]  =  0,  Q  [l]=  1 

Pj:  Ko'Cl]  =  K1’Cl]  =  Q'Cl]=0 

b4:  komCi]-k,*ti]-q-ch-o 
Here  K  '  “  K  "  =  4  etc. 

o  dx  o  cbr 

Evaluation  of  K 
_ Q 

Since  K  [0]=K  Cl3=  RT  '  C  1 3  =  K  "  [t]  =  K  "'[1]  =  0, 

O  O  O  O  O 

K  =  O  X  (X  -  l)4  (10) 

O 

Also  since  K'  [0]=  6,  or  =  6.  Then 

K  =  6x  (x  -  l)4  (11) 

O 

Evaluation  of  Kj 

Since  Kt  [0]  =  [0]  =  Kj  [1]  =  8^  [1]  =  Kj*'  [l]=  0, 
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Kj  =  ex2  <x  -  i)3 


(12) 


Also,  since  Kj"' [  1  ]  =  3,  0  =  y.  Then 

K,  =|  x2  (x  -  l)3 

Evaluation  of  Q 

Since  Q*  CO]=Q’  Cl]  =  Q"  [l]=Q,M  Cl]  =  0, 

Q'  =  Y  x  (x  -  l)3 


and 


Q 


+  C 


From  Q[0]-0t  C  =  0  and  from  Q[  1]  =  1,  y  =  -20. 
Then 


Q  =  1  -  (x  -  l)4  (4x  +  1) 


Permissible  ranges  of  kQ  and  ITt 

Not  all  combinations  of  IT  and  kt  give  desirable  shapes.  It  is 

o 

analyze  possible  limitations  in  terms  of  simple  criteria: 

1.  Zero  condition:  y  =  0. 

Negative  values  of  y  are  meaningless. 

2.  Unity  condition:  y  =  1. 

Bulges  above  y  =  1  are  undesirable. 

dy 

3.  Maximum  or  minimum  condition:  -sr  =  0. 

dx 

Maxima  or  minima  other  than  at  x  =  1  are  undesirable. 

dV 

4.  Inflection  point  condition:  -  =  0. 

dx2 


(13) 


(14) 


(15) 


(16) 


interesting  to 
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Inflection  points  are  considered  undesirable  on  noses  of  bodies. 

As  explained  in  Reference  1,  an  envelope  curve  may  be  determined  for  each 

of  the  preceding  conditions.  As  shown  in  Figure  2,  desirable  values  of  k  and  fT, 

o  1 

are  on  the  inside  of  the  envelope  curve. 

1.  Zero  condition 


y,-hht 


0  ^  x  £  l 


(17) 


The  envelope  in —  and  kj  with  x  as  the  variable  parameter  is  given  by 
k 

o 


h-=!ii=o 

d  X 


(18) 


The  two  envelope  conditions.  Equations  (17)  and  (18)  provide  two  simultaneous 
1  r-' 

equations  In  —  and  kt  which  are  solved  by  the  Cramer  rule  to  give 
K 


1  (x  -  l)4  (3x  +  2)  +  5x  -  2 

/■v 

k 


6x  (x  -  l)4 


(19) 


and 


2(x  -  l)5  -  2(5x  -  1) 

5  = - 

x2(x  -  l)3 


(20) 
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For  x  =  0,  —  =  0  and  =  20. 


For  x  =  1, - ►*>  and  k«  — *  ® 

ir 

o 

The  envelope  curve  is  shown  in  Figure  2.  Desirable  values  of - and  kj  are  on  the 

k 

o 

"inside  curved”  side  of  the  envelope  curve. 

2.  Unity  Condition 


y3  =  h 


c;4r .  c; 


L  kO 


= 1  0^X51 


(21) 


The  envelope  in —  and  Efj  with  x  as  the  variable  parameter  is  given  by 

(T 

o 


(h  -  1)  =  =  0 

dx  *x 


(22) 


The  two  envelope  conditions,  Equations  (21)  and  (22),  provide  two  simultaneous 


equations  in —  and  kt  which  are  solved  by  the  Cramer  rule  to  yield 

(T 

o 


1  3x  +  2 


o 


(23) 
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and 


2(x  -  l)2 


(24) 


For  x  =  0, - and  kx  — >® 

k 

o 

1  5 

For  x  =  1 ,  —  =  —  and  kt  =  0. 
k 

o 

The  envelope  curve  is  shown  in  Figure  2  with  desirable  values  of  —  and  k^on  the 

k 

o 


"inside  curved"  side. 

3.  Maximum  or  Minimum  Condition 


j>y  __ 

dx 


*1 


(25) 


The  envelope  curve  in—  and  fq  with  x  as 


the  variable  parameter  is  given  by 


h"  =  0 


(26) 


The  two  envelope  conditions,  Equations  (25)  and  (26),  provide  two  simultaneous 


equations  in —  and  kx  which  are  solved  by  the  Cramer  rule  to  produce 
k 
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1 


(27, 


5x2 


k  102  -  5x  +  1 
o 

and 

20(x  -  l)2 

- -  (28) 

10x2  -  5x  +  1 


1 

For  x  =  Of  — 
Ic 

o 


1 

Forx  =  lf  — 
k 

o 


The  envelope  curve  is  shown  in  Figure  2  with  desirable  values  of  —  and  E 

IT 

O 

on  the  "inside  curved"  side. 

4.  Inflection  Point  Condition 

For  y3  the  inflection  point  condition  becomes 


3hh"  -  2(h')2  =  0  (29) 


The  envelope  curve  in  —  and  iTj  with  x  as  the  variable  parameter  is  given  by 

o 


=  0,  K",  =  20. 


|.tr1=o. 


3hhf,f  -  h’h”  =  0 


(30) 


A 

These  two  conditions  provide  two  simultaneous  equations  in  —  and  Kt  in  terms  of  x. 

k' 

o 
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X  -w 

Since  the  equations  are  quadratic  in  —  and  k1# 


numerical  procedures  are  necessary. 


The  results  are  shown  in  Figure  2. 

Relative  Fullness  of  Shape 

The  ’’cubic”  polynomial  as  defined  in  normalized  coordinates  specifies  an 
annular  volume  which  is  also  governed  by  the  relative  diameter  of  the  flat  face. 
Consequently  the  prismatic  coefficient  is  a  function  of  not  only  the  ’’cubic”  polyno¬ 
mial  parameters  but  also  the  relative  diameter  of  the  flat  face. 

A  simple  measure  of  the  fullness  of  the  volume  specified  by  the  ’’cubic” 
polynomial  is  given  by  the  integral  F  where 


F  = 


(31) 


For  Equation  (19) 


(32) 


For  constant  values  of  F.  straight  lines  are  indicated  in  Figure  2.  Fuller  shapes 
are  specified  in  the  lower  right  hand  side  of  Figure  2. 


ONE -PARAMETER  ’’CUBIC”  POLYNOMIALS 

A  one-parameter  family  of  ’’cubic”  polynomials  may  also  be  considered  if 
£q9  say,  is  taken  as  the  single  adjustable  parameter.  The  result  is  a  fourth- 
degree  polynomial  or  quartic.  Then  a5  of  Equation  (8)  is  zero  and 

”kj  -  -  12  / ~\  +  8  (33) 

\V 


14 


This  is  plotted  in  Figure  2  as  well  as  the  values  for  the  zero -parameter  "  cubic” 


polynomial,  Equation  (1),  kt  -  2  and — 

k  ^ 


1 

2  * 


By  use  of  Equation  (33)  in  Equation  (9),  the  one-parameter  "cubic"  polynomial 
becomes 


y3=^rK  Cx]+Q[x] 
k  o 
o 


(34) 


where 


K  =  -  6x  (x  -  l)3 
o 


(35) 


and 


Q  =  x2  (3x*  -  8x  +  6) 


A  study  of 


j  permissible  range  of  —  gives 
ko 


(36) 


to  satisfy  the  zero,  unity,  maximum  or  minimum,  and  inflection  point  conditions. 
Representative  curves  are  plotted  in  Figure  3 . 
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